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ABSTRACT 

We investigate the intrinsic opacity of high-redshift galaxies to outgoing ionising 
photons using high-quality photometry of a sample of 27 spectroscopically-identificd 
galaxies of redshift 1.9 < z < 3.5 in the Hubble Deep Field. Our measurement is 
based on maximum-likelihood fitting of model galaxy spectral energy distributions — 
including the effects of intrinsic Lyman- limit absorption and random realizations of in- 
tervening Lyman-series and Lyman-limit absorption — to photometry of galaxies from 
space- and ground-based broad-band images. Our method provides several important 
advantages over the methods used by previous groups, including most importantly 
that two-dimensional sky subtraction of faint-galaxy images is more robust than one- 
dimensional sky subtraction of faint-galaxy spectra. We find at the 3<r statistical con- 
fidence level that on average no more than 4 % of the ionising photons escape galaxies 
of redshift 1.9 < z < 3.5. This result is consistent with observations of low- and 
moderate-redshift galaxies but is in direct contradiction to a recent result based on 
medium-resolution spectroscopy of high-redshift (z w 3) galaxies. Dividing our sam- 
ple in subsamplcs according to luminosity, intrinsic ultraviolet colour, and redshift, we 
find no evidence for selection effects that could explain such discrepancy. Even when 
all systematic effects are included, the data could not realistically accomodate any 
escape fraction value larger than « 15%. 

Key words: galaxies: formation — cosmology: observations — cosmology: diffuse 
radiations 



1 INTRODUCTION 

The diffuse ultraviolet background is one of the key ingre- 
dients in the recipe leading to galaxy formation (Miralda- 
Escude & Ostriker 1990; Giroux & Shapiro 1996). The in- 
tensity of the diffuse ultraviolet background has been esti- 
mated at moderate and high redshifts using the proximity 
effect on QSO absorbers (Liske & Williger 2001; Scott et 
al. 2001; Fernandez-Soto et al. 1995; Kulkarni & Fall 1993), 
and results obtained so far suggest that it is too large to 
be explained by QSOs alone (Madau & Shull 1996, but see 
also e.g. Giallongo et al. 1996). Although systematic errors 
may make these measurements uncertain (Loeb & Eisen- 
stein 1995, Pascarelle et al. 2001), several suggestions have 
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been put forward to explain the observed extra intensity. 
One of the most popular suggestions calls for the escape of 
a large fraction of ionising photons from high-redshift galax- 
ies. If correct, this would imply that the interstellar media of 
high-redshift galaxies must be at least partially transparent 
to ionising photons. 

But the best measurements available at low and moder- 
ate redshifts show that galaxies are (at least nearly) opaque 
to ionising photons. The large column densities of neutral 
hydrogen that surround most galaxies should be enough, 
depending on how they are distributed, to very effectively 
quench any outgoing flux of ionising photons. In fact, Lei- 
therer et al. (1995) find, in a sample of four starburst galax- 
ies at an average redshift (z) ~ 0.02, that the limits to 
the escape fraction of Lyman limit photons range from 
/esc < 0.0095 to / csc < 0.15 (with an observed flux ratio 
at 900 A and 1500 A of F 9/15 < 0.10 - 0.20). Hurwitz et al. 
(1997) analyse the same data and obtain slightly less strin- 
gent limits, ranging from / csc < 0.032 to / csc < 0.57. Dehar- 
veng et al. (2001) find a more stringent limit / osc < 0.064 
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for a starburst galaxy of redshift z = 0.0448. At higher red- 
shifts [z fa 1), Ferguson (2001) finds similar upper limits 
/esc < 0.20 to the escape fraction. 

It is much less clear whether galaxies at still higher 
redshifts are as opaque to ionising photons as galaxies at 
low and moderate redshifts. Accurate measurements be- 
come progressively more difficult at higher redshifts, because 
galaxies are much fainter and the presence of Lyman a forest 
absorption introduces further uncertainties. The only two 
measurements available so far are from Steidel, Pettini & 
Adelberger (2001) and Giallongo et al. (2002). Adopting a 
correction factor for intervening Lyman a absorption deter- 
mined from a composite QSO spectrum, Steidel et al. re- 
ported F g/15 = 0.22 ± 0.05 using 29 galaxies at (z) = 3.40, 
while Giallongo et al. reported a lcr upper limit -Fg/15 ^ 0.05 
using two galaxies at z ~ 3. The former measurement indi- 
cates that galaxies at 2 ~ 3 are much more transparant to 
ionising photons, contributing an equal amount of ionising 
flux to the ultraviolet background radiation as the QSOs. 
The latter measurement indicates otherwise. 

In this article, we present a new measurement of / csc 
using 27 galaxies at redshifts 1.9 < z < 3.5 in the Hubble 
Deep Field (HDF). All of these galaxies have secure spec- 
troscopic redshifts (see Cohen et al. 2000-COO-and Dawson 
et al. 2001-D01) and accurate broad-band photometry from 
HST/WFPC2 observations in the F300W, F450W, F606W, 
and F814W bandpasses. Our measurement is based on deep 
space-based images, which provides two important advan- 
tages over the methods used by previous groups: First, our 
sky subtraction is more accurate, because two-dimensional 
sky subtraction in faint-galaxy images is more robust than 
one-dimensional sky subtraction in faint-galaxy spectra. Sec- 
ond, the wavelength interval below the rest-frame Lyman 
limit over which we integrate is larger, because the space- 
based broad-band images are sensitive to observed-frame 
wavelengths as short as A ~ 2800 A, allowing us to inte- 
grate over rest-frame wavelength intervals as large as « 100 
A in some cases. Our method also improves upon previous 
work by accounting for the effect of the Lyman a forest ab- 
sorption by performing a large numer of random realizations 
of a parameterized distribution, to account for variations of 
the absorbing clouds across different lines of sight. 

Our results show that the escape fraction estimated 
from the ensemble of 27 galaxies has a (statistical) 3a up- 
per limit of /esc < 0.04. Given that our galaxy sample spans 
a range of properties in redshift, intrinsic colours, and lu- 
minosities, we also explore possible systematic effects by 
dividing the galaxies into different subsamples. We do not 
find any evidence to support that bluer galaxies have higher 
escape fractions of ionizing photons. Instead, we find that 
more luminous galaxies appear to have shallower ultraviolet 
spectral slopes and higher escape fractions (at the 2a level 
of significance) . We do also search for other possible system- 
atic effects that could be inherent to our method, and find 
that other effects cannot in any realistic way account for a 
large increase in the escape fraction. 

2 DATA 

We consider all galaxies in the HDF with known spectro- 
scopic redshifts in the range 1.9 < z < 3.5. The lower limit 



Table 1. Properties of the galaxy sample used in the analysis. 
Celestial coordinates arc J2000. 
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is chosen so that the rest-frame Lyman limit begins to af- 
fect the photometry in the F300W bandpass. The upper 
limit is chosen so that we can still reasonably distinguish 
absorption due to the Lyman a forest from the intrinsic Ly- 
man limit absorption. We previously measured high-quality 
photometry of these galaxies using HST/WFPC2 through 
the F300W, F450W, F606W, and F814W filters (Fernandez- 
Soto, Lanzetta & Yahil 1999-FLY99 hereafter; Lanzetta et 
al. in preparation) . The sample includes 27 galaxies, as given 
in Table 1. The reader is referred to FLY99, COO, D01, and 
Fernandez-Soto et al. 2001 for more details on individual 
galaxies. 

We determine the absolute AB magnitude of each 
galaxy using a fiducial spectral energy distribution deter- 
mined from a template fitting technique that we will describe 
in the next section. We adopt a cosmological model with vac- 
uum energy density Qa = 0.65, matter density Q.m = 0.35, 
and a Hubble constant Ho = 65 km s _1 Mpc -1 . Absolute 
AB magnitudes of the galaxies at rest-frame wavelengths 
A = 1500 A range from w -21.8 to ^ -25.4. Apart from 
these absolute magnitudes, all other results in this paper are 
independent of the chosen cosmological model. 



3 METHOD 

Our method is based on maximum-likelihood fitting of 
model galaxy spectral energy distributions, including the 
effects of intrinsic Lyman-limit absorption and random re- 
alizations of intervening Lyman-series and Lyman-limit ab- 
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sorption, to space- and ground-based broad-band images. In 
practice, our method proceeds in several steps as follows: 

First, we determine a fiducial spectral energy distribu- 
tion of each galaxy by fitting model spectral energy distri- 
butions to the broad-band photometric measurements that 
fall between rest-frame wavelengths A = 1250 and 2800 A. 
The wavelength range is chosen to avoid the spectral dis- 
continuity produced by the Lyman series (rest-frame wave- 
length A < 1215 A) and the Balmer break (rest-frame wave- 
lengths A 3000 A). The fiducial spectral energy distri- 
butions are based on the six spectrophotometric templates 
(of E, Sbc, Scd, Irr, SB1, and SB2 galaxies) of our previous 
photometric redshift measurements (Yahata et al. 2000; see 
also Bem'tez 2000) together with a power-law form f v cx \ a . 
(The power-law form is required because a few of the galax- 
ies are found to be bluer-in the ultraviolet rest-frame-than 
any of the starburst spectrophotometric templates. In each 
of these cases, a < 0.5.) The spectral energy distributions 
are selected according to a \ 2 analysis similar to the analy- 
sis used for our previous photometric redshift measurements 
(Lanzetta et al. 1996; FLY99). 

Next, we modify the fiducial spectral energy distribu- 
tion of each galaxy by adding the effect of intrinsic Lyman- 
limit absorption, which is characterized by an assumed 
Lyman-limit optical depth r g that is taken to be univer- 
sal across all galaxies. For a Lyman-limit optical depth r g , 
the output spectrum /(A, r g ) in terms of the input spectrum 
/(A, 0) at rest-frame wavelength A is 

/(A, Tg ) = /(A, 0) exp[-T g (A/A LL ) 3 ] (1) 

at rest-frame wavelengths A < All, where All is the rest- 
frame wavelength of the Lyman limit. 

Next, we modify the fiducial spectral energy distribu- 
tion of each galaxy by adding the effect of random realiza- 
tions of intervening Lyman series and Lyman limit absorp- 
tion. Here we characterize the distribution of Lyman a for- 
est absorption systems by a distribution function F(N, b, z), 
which is defined in such a way that F(N, b, z) dN db dz is 
the expected number of absorption systems in the neutral 
hydrogen column density interval dN around N, Doppler pa- 
rameter interval db around b, and redshift interval dz around 
z. We adopt as the functional form of the distribution func- 
tion 

F(N,b,z) cc {l + zyN-^G*{b,a b ,b miYl ), (2) 

where G* (b , ) is a "truncated" Gaussian distribution 

of mean b, dispersion at, and minimum truncation value 
bmin- We adopt 7 = 2.5, /3 = 1.6, b — 23 km s _1 , a b = 8 km 
s~\ and & m i n = 15 km s _1 . 

We have extensively checked that simulated spectra 
generated according to these parameters reproduce to within 
observed scatter the average Lyman a absorption Da and 
average Lyman series absorption Db- In Figure 1 we show 
the average values of Da and Db, compared to values found 
in the literature. The agreement is perfect within the scatter 
present in the observational data. 

The incidence of Lyman limit systems (those with hy- 
drogen opacity r > 1) is an important factor in our simula- 
tion, as they are amongst the main drivers of the absorption 
bluewards of the galactic Lyman edge. We plot in Figure 2 
the density of Lyman limit systems generated by our Lyman 
a forest model, and compare it with the numbers measured 
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Figure 1. Values of the Lyman decrements Da and Db obtained 
by averaging 1000 spectra generated by our model at each red- 
shift (continuous line). The curves are compared with the values 
by Webb (unpublished, dotted line) which our group uses for pho- 
tometric redshift estimation, and with a compilation of Da values 
obtained from the literature (Okc & Korycansky 1982; Bcchtold et 
al. 1984; Stcidel & Sargent 1987; O'Brian, Wilson & Gondhalekar 
1988; Schneider, Schmidt, & Gunn 1989; Giallongo & Cristiani 
1990; Schneider, Schmidt, & Gunn 1991). 



by different authors (as described in the Figure). As can be 
seen, our model produces slightly less Lyman limit systems 
at the lowest redshifts than have been observed. Whereas 
this effect could produce a bias favouring lower escape frac- 
tion values, we have checked that this effect is not important 
(see SS4.2 and 4.3 for details). This slight deflect is the re- 
sult of our decission to adopt a single parameterization to 
generate all of the Lyman a forest-a change in the normali- 
sation to improve the aspect of Figure 2 would conflict with 
the data shown in Figure 1. Globally, we consider that our 
model reproduces within reasonable limits the properties of 
the intergalactic medium. 

Next, we measure the simulated spectrum of each 
galaxy — constructed from the fiducial spectral energy dis- 
tribution of each galaxy modified by the effects of intrin- 
sic Lyman-limit absorption with an assumed Lyman-limit 
optical depth r g and of random realizations of interven- 
ing Lyman series and Lyman limit absorption — at short- 
wavelengths, i.e. over the F300W bandpass which is sensi- 
tive at wavelengths below the rest-frame Lyman limit, and 
we compare observed energy fluxes f^j g and uncertainties 
^obs an d simulated energy fluxes of the ith galaxy in 

the jth simulation. To randomize the effects of intervening 
Lyman series and Lyman limit absorption, we repeat this 
procedure a large number of times, and we form the likeli- 
hood as a function of Lyman-limit optical depth r g as 
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Figure 2. Number density of Lyman limit systems as observed 
by several authors (see labels) and produced by our method. The 
difference between our model and the average of the different 
observations is always less than 0.2dex in the redshift range of 
interest, and almost zero at z > 2.8. 
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where the product extends over the JV ga i galaxies of the sam- 
ple and where the sum extends over the N s i m simulations of 
the analysis. We find that AT Bim = 100 is enough to represent 
the variation of intervening Lyman series and Lyman limit 
absorption. 



4 RESULTS AND DISCUSSION 

Figure 3 shows results of the method described in § 3 applied 
to five examples of the galaxies. The examples are chosen 
to span a range in spectral type, redshift, and luminosity. 
For clarity, only 10 of the 100 realizations are presented for 
each galaxy, ft is obvious from inspection of Figure 3 that 
in many cases (especially at the lowest redshifts) the fluxes 
observed through the F300W filter are less than the fluxes 
expected, if the galaxies were completely transparent to ion- 
ising photons. This trend is also present in the other galaxies 
of the sample. Our quantitative maximum-likelihood analy- 
sis of this trend across the entire sample of galaxies estab- 
lishes a lower limit to the optical depth at the Lyman limit 
or an upper limit to the escape fraction for these galaxies. 
In this section, we discuss these limits. 



determined from the entire sample of galaxies. The best-fit 
value of the galaxy Lyman-limit optical depth determined 
from the entire sample of galaxies is 



5.2^ 



(4) 



and the corresponding best-fit value of the galaxy escape 
fraction determined from the entire sample of galaxies is 



0.008 ± 0.006. 



(5) 



These results apply at the average redshift of the entire sam- 
ple of galaxies, which is 



(z) = 2.77 



(6) 



The 3<r upper limit to the galaxy escape fraction determined 
from the entire sample of galaxies is 



/ csc < 0.039. 



(7) 



Apparently, our analysis rules out the possibility of a large 
galaxy escape fraction at high significance but cannot rule 
out (at the 3a level) the possibility of a galaxy escape frac- 
tion identically equal to zero. Based on these results, we 
conclude that galaxies of redshift z « 3 are highly opaque to 
ionising photons. 



4.2 Possible selection Effects 

Results of our analysis are apparently at odds with results of 
Steidel et al. (2001) that galaxies of {z) « 3.4 are not highly 
opaque to ionizing photons. These authors observed in the 
average spectrum of 29 galaxies a flux ratio F g / 15 = 0.056. 
After correcting for the effects of the Lyman a forest absorp- 
tion using a composite QSO spectrum at the same redshift, 
they found an intrinsic flux ratio of F g / 15 — 0.22 ± 0.05- 
where the effects of intrinsic galactic absorption and attenu- 
ation by dust or gas have not been eliminated. The authors 
further incorporated to this measurement a plausible esti- 
mate of the flux discontinuity at the Lyman limit due to 
absorption in stellar atmospheres to derive an escape frac- 
tion /esc^O.10. 1 

Steidel and collaborators suggested that the high escape 
fraction found in their work might result from some selec- 
tion bias of their galaxy sample. The galaxies incorporated 
into their analysis include only 29 of 875 galaxies for which 
they obtained spectroscopic observations with the Keck tele- 
scope, selected on the basis of quality of the observations at 
blue wavelengths. This selection condition quite naturally 
selects galaxies that are either very blue, very luminous, or 
both. QSOs are known to be transparent to ionizing photons 
at ultraviolet wavelengths (Zheng et al. 1997), so it is not 
unreasonable to speculate that very blue or very luminous 
(or both) galaxies might exhibit larger escape fractions than 
avarage galaxies. To examine this possibilty, we repeated the 
analysis described in § 3 for various subsamples of the galaxy 
sample. 



4.1 Galaxy Escape Fraction 

Figure 4 and Table 2 show results of the method described 
in § 3 applied to the data described in § 2. The first panel of 
Figure 4 shows the relative likelihood as a function of log(r g ) 



Steidel and collaborators used a definition of / eS c which is dif- 
ferent from the usual one— they normalised the escape fraction at 
900 A to the escape fraction at 1500 A. They reported / esc 0.50, 
which is converted here to the usual definition assuming the es- 
cape fraction at 1500 A is ss 0.15 — 0.20 as quoted by the authors. 
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Figure 3. Photometry and model spectra (10 simulations each time) of a sample of galaxies with no intrinsic Lyman limit absorption 
(r g = 0, left) and strong intrinsic absorption (r g = 10, right). The vertical lines mark the positions of the Lyman limit, Lyman (3 and 
Lyman a lines, and A = 2800A in the rest frame. Circles with error bars show the photometric measurements. The squares show the 
photometry obtained from the fitted spectrum+HI absorption. The galaxies have been chosen to show the full range in redshift, type, 
and luminosity. 



Color: First, we repeated the analysis for two subsam- 
ples of the galaxy sample divided according to color. Re- 
sults are shown in the second panel of Figure 4 and the 
second group of entries in Table 2. Here "bluer" galaxies are 
galaxies best described by SB1 or power-law spectrophoto- 
metry templates and "redder" galaxies are galaxies best de- 
scribed by Scd, Irr, and SB2 spectrophotometric templates. 
Results obtained from both subsamples are consistent with 
each other and with results obtained from the entire sample 
of galaxies, although there is a (statistically insignificant) 
trend for the bluer galaxies to indicate a smaller galaxy es- 
cape fraction than is indicated by the redder galaxies. This 
result runs contrary to the idea described above. 

Luminosity: Next, we repeated the analysis for two sub- 
samples of the galaxy sample divided according to luminos- 
ity. Results are shown in the third panel of Figure 4 and 
the third group of entries in Table 2. Here "low-luminosity" 
galaxies are galaxies of absolute magnitude AB > —23.75 
and "high-luminosity" galaxies are galaxies of absolute mag- 
nitude AB < —23.75. Results obtained from both subsam- 
ples are consistent with each other and with results obtained 
from the entire sample of galaxies, although there is a (statis- 



tically insignificant) trend for the lower-luminosity galaxies 
to indicate a smaller galaxy escape fraction than is indicated 
by the higher-luminosity galaxies. 

Redshift: Finally, we repeated the analysis for two sub- 
samples of the galaxy sample divided according to redshift. 
Results are shown in the third panel of Figure 4 and the third 
group of entries in Table 2. Here "low-redshift" galaxies are 
galaxies of redshift 1.95 < z < 2.85 and "high-redshift" 
galaxies are galaxies of redshift 2.85 < z < 3.50. Results ob- 
tained from both subsamples are consistent with each other 
and with results obtained from the entire sample of galaxies, 
although there is a (statistically insignificant) trend for the 
lower-redshift galaxies to indicate a smaller galaxy escape 
fraction than is indicated by the higher-redshift galaxies. To 
some extent, this statistically insignificant difference could 
be due to the already mentioned relative lack of Lyman limit 
systems in our model at low redshift. Given that the results 
from the high-redshift sample (where our model reproduces 
exactly the observed number of Lyman limit systems) are 
not significatively different from the results obtained using 
the low redshift sample, we infer that the difference in the 
model does not induce a large change. 
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Figure 4. Relative likelihoods of the value of r g for the complete 
sample (left panel). The other panels show the results for differ- 
ent subsamplcs, broken according to (left to right) type, rcdshift, 
and luminosity. In all panels the horizontal lines correspond to 
statistical 1, 2, and 3 a confidence intervals. 



It should be remarked that none of the subsamples anal- 
ysed is different from the values given by the complete sample 
at the 3a level. Neither are any of the subsamples incompati- 
ble with their complementary subsamples at the same level. 
With this in mind, it is nevertheless interesting that the 
apparent effects of luminosity and ultraviolet blueness on 
/ csc seem to pull in opposite directions, suggesting that the 
higher escape fraction reported by Steidel and collaborators 
is unlikely due to selection effects. 



4.3 Systematic uncertainties 

All the confidence limits reported above refer only to the sta- 
tistical properties of the sample of galaxies we are analysing. 
It is very likely that the method we are using to analyse the 
data may introduce errors of a systematic nature, that could 
potentially be of the same order of the statistical ones. We 
have performed some checks in order to estimate which could 
be the contribution from several sources of uncertainty. 



4-3. 1 Uncertainties induced by the Lyman a forest model 

As was shown in §3 above, our Lyman a forest model has 
been checked against the observed properties, and agrees 
with them within observational limits. In order to estimate 
the uncertainty in / csc that could originate from the uncer- 
tainties in the model, we have performed new measurements 
of /esc using forest models that differ from our "standard" 
one by applying a 10% increase/decrease in the normali- 
sation of the line density. We must remark that this 10% 
change does largely overestimate the possible uncertainties 
in our model-the models with 10% more or less lines repre- 
sent bad fits to the observed forest properties. 

As expected, an increase in the forest density implies an 
increase in the measured value of the escape fraction: more 
of the UV photons can escape each galaxy while keeping the 
observed photometric properties, as their mean free path 
is reduced by the increased density of absorbers. Mutatis 
mutandis, when the forest density is reduced, the value of 



/ osc moves down. The change is, however, small: at the 3a 
confidence level, the original limit (/ esc < 0.039) becomes 
/csc < 0.030 (/esc < 0.046) when the line density is decreased 
(increased) by 10%. 

Another potential problem in our forest model could be 
the lack of clustering in our simulations, as opposed to what 
is shown by the observations. Our measurement is not af- 
fected by the small-scale clustering that has been observed 
in the Lyman a forest lines (Fernandez-Soto et al. 1996) 
because we are integrating over much larger scales. Of po- 
tential importance is the observed clustering of high-column 
density absorbers (Sargent, Boksenberg & Steidel 1988). It 
is difficult to quantify its effect, but we expect it should be 
dilluted by the size of our sample, and estimate that it pro- 
duces a smaller variance than the one induced by the ±10% 
change in the number of lines presented above. 



4-3.2 Uncertainties induced by the choice of the spectral 
templates 

In the case of some of the galaxies, the photometric data 
available could allow for the assignation of two different 
fiducial spectral templates with similar goodness-of-fit val- 
ues. The choice of one of them over the other (which we 
perform algorithmically via the minimum chi-squared cri- 
terium) could produce an extra uncertainty in the measure- 
ment of /esc • The same effect can be produced if our selection 
of templates is not dense enough to cover the possible UV 
shapes, or if the photometric uncertainties are too big. 

We have evaluated the uncertainty induced by this ef- 
fect by performing the following exercise: we have forced all 
spectral templates to be one step bluer (always along the se- 
quence formed by Ell, Sbc, Scd, Irr, SB2, SB1) than the one 
which is the best-fit choice 3 , and repeated the analysis. This 
does in fact result in very bad fits for many of the galaxies, 
and must thus be considered a very large overestimate of the 
possible uncertainties. Under this scenario, with the galaxies 
emitting far more ionising photons than what is indicated 
by their actual UV continua, the escape fraction diminishes 
to very low values: / CS c < 0.016 at the 3a level. 

We have also performed the opposite exercise, forcing all 
galaxies to be one step redder than their best-fit template 4 . 
Under this assumption the galaxies emit far less ionising 
photons than what their UV continua actually suggests, and 
this results in an increase in the acceptable escape fraction, 
which reaches / OS c = 0.065 , with a 3a upper limit / CS c < 
0.156. 

As we explained before, these limits must be taken as 
very large overestimates of the possible systematic effects 
induced by our spectroscopic templates, as we do not realis- 
tically expect to mistake all galaxies in the same direction. 
However, even under this extreme assumption, and incorpo- 
rating also the effects described in §4.3.1, the upper limit 
to the measured escape fraction is still far below the values 
measured by Steidel et al. (2002). 



3 Exception made of the ones that are already fitted by a SB1 
template or a power-law 

4 Exception made of those galaxies fitted by a power-law 
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Table 2. Measurements of / esc obtained via the maximum likelihood method described in the text. The different subsamples used in 
the analysis are listed, together with the results and statistical confidence intervals associated to each one. 



Sample 


N 2 


/esc 


la Interval 


2<t Interval 


3<t Interval 


All 


27 


0.008 


(0.001-0.013) 


(0.000-0.023) 


(0.000-0.039) 


Redder 


15 


0.008 


(0.004-0.020) 


(0.000-0.036) 


(0.000-0.054) 


Bluer 


12 


0.000 


(0.000-0.003) 


(0.000-0.018) 


(0.000-0.047) 


High Luminosity 


13 


0.026 


(0.017-0.050) 


(0.008-0.072) 


(0.002-0.093) 


Low luminosity 


11 


0.000 


(0.000-0.001) 


(0.000-0.004) 


(0.000-0.013) 


Low redshift 


13 


0.004 


(0.000-0.008) 


(0.000-0.018) 


(0.000-0.030) 


High redshift 


11 


0.169 


(0.056-0.389) 


(0.008-0.857) 


(0.000-1.000) 



5 CONCLUSIONS 

We have presented a new method to measure the escape 
fraction of ionising photons from galaxies at z ~ 3. Our 
method provides a cleaner measurement compared to the 
classical spectroscopic technique because the problem of the 
sky subtraction is reduced almost to insignificance and we 
can use a wider bandpass. 

The value we obtain (/ csc = 0.008, with / csc < 0.039 at 
a 3<r confidence level) is in agreement with most previous 
measurements of the escape fraction of ionising photons at 
high and low redshift, being more stringent than most previ- 
ous upper limits. It is largely in contradiction with the value 
presented by Steidel et al. (2001) for high-redshift galaxies . 
We have studied possible selection effects that could cause 
this apparent difference and conclude that, within the limits 
set by the size of our sample, they cannot explain the large 
discordance between our results. Neither can the difference 
be explained by possible systematic effects induced by our 
method, as we also have shown. A (very conservative) es- 
timate of the confidence limits when all systematic effects 
are included leads to a 3a upper limit to the escape fraction 
that cannot be higher than / eac 0.15. 

If the escape fraction of ionising photons reported here 
represents the general situation at high redshift, then normal 
galaxies cannot be responsible for any significant fraction of 
the high-redshift ionising background. On the other hand, 
if the proximity effect measurements of the background flux 
are confirmed, then the problem persists to find the objects 
responsible for the ionisation state of the high-redshift uni- 
verse. Deep, narrow-band imaging of local and high-redshift 
galaxies at wavelengths slightly above and below their intrin- 
sic Lyman limits could settle this argument. Such imaging 
campaigns would be less expensive, in terms of observing 
time, than the spectroscopic observations done to date. 
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